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Abstract 

We study the bulk viscosity of a pion gas in unitarized Chiral Perturbation Theory at low and moderate temperatures, 
below any phase transition to a quark-gluon plasma phase. 

We argue that inelastic processes are irrelevant and exponentially suppressed at low temperatures. Since the system 
falls out of chemical equilibrium upon expansion, a pion chemical potential must be introduced, so we extend the 
existing theory to include it. We control the zero modes of the collision operator and Landau's conditions of fit when 
solving the Boltzmann equation with the elastic collision kernel. 

The dependence of the bulk viscosity with temperature is reminiscent of the findings of Fernandez- Fraile and Gomez 
Nicola [1] , while the numerical value is closer to that of Davesne [2] . In the zero-temperature limit we correctly recover 
the vanishing viscosity associated to a non-relativistic monoatomic gas. 
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1. Motivation 

Bulk viscosity C is responsible for the equilibration of a system subject to dilatation or compression per- 
turbations. Although it usually is much smaller than the more common shear viscosity rj, it is an interesting 
probe of the loss of dilatation invariance [3] through the trace anomaly in quantum field theory. In the 
context of relativistic heavy ion collisions, bulk viscosity can be experimentally accessed for example by 
means of fluctuations of the two-point correlation of the energy-momentum tensor [4] . It has been computed 
in the asymptotically high temperature regime of Quantum Chromodynamics [5] , and found to be governed 
there by inelastic processes. 

However, published investigations concerning bulk-viscosity in the low-energy pion gas (the opposite, low- 
temperature limit of strongly interacting matter at zero baryon number) yield contradictory results. While 
the Chiral Perturbation Theory based computation of Fernandez-Praile and Gomez Nicola features a clear 
maximum of C associated with the scale-violating pion mass [1], this is absent from earlier studies [2,6] based 
on empirical phase-shifts, that also differ in the numeric values of C- 

Neither of these calculations carefully identifies the zero-modes of the collision operator that might change 
the temperature counting of the viscosity. A chiral limit calculation has however been reported [7] that studies 
these zero modes in elastic scattering. 

Further, none of these works properly emphasizes the role of the Landau-Lifschitz constraint, that is 
however properly treated in a linear sigma model computation [8] . 

Given the situation, we feel it is worth to settle the issue of the bulk viscosity of a low-temperature pion 
gas by adopting the best features of the various existent approaches. In this article we provide a Unitarized 
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ChPT-based computation at the physical pion mass m^, where the interaction kernel reproduces the physical 

phase-shifts, and where the zero modes and Landau-Lifschitz constraints are correctly implemented. 

Since we find that elastic processes dominate the low-temperature bulk viscosity, pion number is conserved 
in the relevant time-scales, and the need to introduce in the theory a pion chemical potential /i arises, which 
is our additional contribution to the discussion. We also want to emphasize that the rclativistic chemical 
potential satisfies fi < m,r instead of /x < in the non-relativistic theory. Thus, some of prior work with 
H = might actually not be easily mapped to the non-relativistic limit (where it would correspond to an 
artificial /i = — m^r if care was not exercised in subtracting ni^ from the energy), finding a divergence of 
CiT) at vanishing temperature that might be purely academic. 

We make no attempt to extend the validity of the approach to higher temperatures where strangeness is 
an active degree of freedom and our considerations shoiild be imderstood as valid up to temperaturcis in the 
range 120-140 MeV at most, although we may plot results at slightly larger temperatures for archival and 
comparison purposes. 



2. Dominance of elastic processes 



Departures from equilibrium of the one-particle distribution function fp = np + Sfp that are associated 
with system expansion can be parametrized, in linear order of velocity gradients, as 

Sfp = -ripil + np)p W-Y Ap, (1) 

with Up the Bosc-Einstcin equilibriimi distribution fimction np = [e^^ — , E the rclativistic kinetic energy, 
and /3 = l/T the inverse temperature. Ap must be an adimensional function of the particle momentum. 
Therefore, defining the dimensionless variable x = EpjmT^ = ^/T+Jp/Tn^, one can expand it as 

Ap~ao + aiEp + fia^^p + • • • (2) 

We can conveniently also consider Ap as a function A(x). From dissipative hydrodynamics, the bulk viscosity 
is known to be a simple integral of Ap 

with g = 3 being the pion isospin degeneracy. Relaxation of fp towards equilibrium in a dilute gas is governed 
by the Boltzmann equation Pfid^^f = C[f] that is linearized in terms of the separation from equilibrium 5f, 
yielding a linear system in an appropriate polynomial basis, for example that of Eq. (2), resulting in C\A) = 
\S). 

The collision operator C has parts associated with elastic tttt — > tttt and inelastic tttt — ^ TTTTTTTr processes, 
C = C2 + Ci + ... 

Writing explicitly the elastic part C2 and assuming for the time being that the system does not have a well 
defined particle number (no chemical potential), the linearized Boltzmann's equation with a compressive 
disturbance has the well-known form 

\s) ^ (1/ - Ey^ = 

7, — 77T — V I dTi2^:ip{l + nx){l + n2)n:inp[Ap + Ap^ Ap^ Ap2 

] + C,\A) (4) 

ZUp^L -T Tip) J 

with (iri2,3p carrying the elastic interaction and defined below in Eq. (17), and v"^ being the squared sound 
velocity dP/dp. 

The zero modes of the elastic collision operator C2 can be identified from the structure Ap+Ap^ —Ap-^—Ap^ 
and correspond to Aq = 1 (associated to particle-number conservation since the number of terms in this 
equation is two positive and two negative due to the 2^-2 scattering) and Aq = x = Ep/m.,^ associated 
to energy conservation (this mode is also present in the inelastic parts of the linearized collision operator). 
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Thus the collision operator must be inverted in a linear space perpendicular to these modes. Assuming 
this has been achieved, the linearized Boltzmann equation returns a solution Apj_. However, this solution is 
not unique as A± + dix will also solve it for arbitrary ai including both collision operators C2 + C4, and 
A± + ao for arbitrary Uq including only C2. These two constants can be specified by the two Landau-Lifschitz 
constraints ( "conditions of fit" ) necessary to fully define the reference frame and the particle density, that 
are 



Substituting dfp in terms of A, these can also be written as 



/(TV f drv 

^^np{l + np)ApEl = ; j j^^M^ + np)ApEp = 



(6) 



The Boltzmann equation can then be symbolically written as 

{C2 + Ci)\A^ + Ao) = \S) or C2\As_) + dHAi,) + \Ao)) = \S) (7) 

in view that ^0 contains the zero modes of C2. 

Now wc seek to see which of these terms dominates the inversion of the collision operator at low tem- 
peratures. It is clear that, for massive particles, one can use non-relativistic power counting for low-enough 
temperature, so that E ~ +p^/2m.„, with m,r oc 0(1), {Ep — m^) oc 0{kBT), p oc \fm^. Then it is 
possible that A\^ be suppressed by a temperature power respect to j4o given that the zero modes involve 
the lowest powers of the energy 1, x. This is however overcompensated by the exponential suppression of 
inelastic processes C^jC^ oc e"^™"/^ from the additional Up factors, and an additional two powers of T from 
chiral suppression of the inelastic amplitude. Therefore, at low temperatures C4IA0) is negligible respect to 
C2\A^. In the following we neglect the contribution of C4, that however should become more and more 
relevant as one approaches the phase transition (or alternatively, when typical CM collision energies are of 
the order of the inelastic resonance /2(1270)). However the strangeness degree of freedom is active before 
that energy, so including C4 without attending to the kaons in the bath is not a sensible approach. 

We therefore concentrate on elastic processes and proceed to analyze C2 in detail. Since elastic collisions 
do not change pion number, there is an approximate conservation law that requires the introduction of a 
chemical potential. 

3. Bulk viscosity in the presence of a pion chemical potential 

In the presence of a pion chemical potential, the left hand side of the Boltzmann equation (4) needs to be 
modified 

p^5^/o(x)|^ = /3np(l + Up) _ ^2^2 _ ^^^-1^ y . ^ 
with the inverse compressibility at fixed density and the sound velocity at constant particle density defined 



as 



where n = n{T, jj) is the particle number density and p = p{T, p) is the energy density. The source function 
in the left hand side of Eq. (4) is therefore to be substituted by (^p^ — Epv'^-^ — EpK'^). 

It is then convenient to employ the Landau-Lifschitz conditions (6) to add two (vanishing) terms to Eq. (3) 
to make the bulk viscosity take an elegant form 

^ = f / j^Ep^^^^ + "^^^^ (y - - ■ 
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Since the collision operator is a linear operator acting on A, it is also convenient to define a scalar product 
in the function space given by the measure 



^^^^^^iJ^^"'^^^'^^^''' ^^^^ 

and a basis of the function space. After much trial and error we have found practical to employ a basis that 
is not totally orthogonal, but where the first and third polynomials of a; = -B/m^ are given by 

W^i; P.(.).(i-.^).^-^.-i. (12) 

All other polynomials Pi (a;), Pzix) and higher are then chosen by standard orthogonalization to these two, 
and monic, that is, in P„ the coefficient of x" (maximum power of x) is 1. The advantages of these basis 
are that the first two polynomials Pq and Pi span the zero modes of the collision operator, and the second, 
P2, is the source 15) function, bringing about certain simplification. (We should like to stress again that for 
i> 1 the P, form an orthogonal, monic basis). 

Equation (10) for the bulk viscosity is then the projection of the solution to Boltzmann's equation over 
the source, 

C = ^{Aix)\P,{x)) . (13) 



4. Solution of Boltzmann equation and satisfaction of constraints 



We seek the coefficients of the expansion 

oc 

A{x) = Y,anPn{x) (14) 

n=0 

that satisfy both the quantum Boltzmann equation C2\A) = \S) and the Landau-Lifschitz constraints 
{E,\A)=0,{E',\A)=0 

By truncating the polynomial basis {Pi(x)} to a finite number of polynomials, one can represent the 
linearized Boltzmann equation as a symmetric linear system that can be inverted on a computer, in terms 
of Cij = {Pi\C\Pj) = Cji that is the matrix element of the collision operator in this basis. 
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This system is compatible but indeterminate. Inclusion of Po{x) = 1 in the space where C is to be inverted 
makes the system incompatible, since it is a zero mode and the source has non-zero projection {Pq{x)\S) ^ 
0. It is not possible to find a basis in which Po(x) and Pi(a;) (spanned by the two zero modes) are both 
orthogonal to P2{x) = \S). In the function space Pi . . . P„ that we keep, the source, by construction, has 
only a projection over P2, S2 = {S\P2)- 

As for the indeterminacy, note that ai is so far completely arbitrary. 

The first step is thus to invert C in a subspace {P2 . . . P„}. If one truncates at n = 2, for example, then 
02 = 02282- 

We then identify a„ = a„, for n > 2, and we fix ao and oi by demanding satisfaction of Eq. (6). For 
example, truncating at n = 2, we would find 
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ao{x^\Po) + ai{x^\Pi) = -a2{x^\P2) (16) 
ao{x\Po) + ai{x\Pi) = -a2{x\P2) 

For higher n it is a matter of solving a more complex 2x2 linear system, task that we perform easily with a 
computer program. Essentially all the computer time is spent in obtaining the Cij matrix elements by means 
of Montecarlo integration over the multidimensional space of particle momenta, with the linear algebra in 
this section costing negligibly small CPU effort. 

Note that the space where we actually solve the Boltzmann equation is not perpendicular to the zero modes 
(our discussion in section 2 above was too naive). However, the inversion is well conditioned and convergence 
in n very fast because there is a non-zero projection of both zero modes that is perfectly perpendicular to 
the inversion function space, so that no eigenvalue ever vanishes. 

The resulting A{x) can then be substituted in either the Boltzmann equation or the Landau-Lifschitz 
conditions and satisfies both uniquely. 



5. Low-energy pion scattering 



The interaction between pions (1 + 2 — )• 3 +p) appears in the squared matrix element |T|^ 

3 

rfrr.,3. = (^(2.)^.(fc, + k2-k,-p) (17) 

averaged over isospin and partial waves. We employ the Inverse Amplitude Method (lAM) [11] to represent 
|Tp. The interaction is exactly the same as in [10] and well documented in the literature, in particular we 
employ the fits from Ref. [12]. The key of the method is to employ a dispersive analysis to describe pion-pion 
scattering in which the fimction in the dispersive integral is not the partial wave itself but its inverse. The 
subtraction constants and left cut are approximated in Chiral Perturbation Theory, while the right cut is 
known exactly in the elastic region (and to a very good approximation up to about 1.2 GeV). The resulting 
lAM formula is then written down in terms of the Chiral Perturbation Theory expansion for the partial 
waves, t ~ ^2 + ^4 + ■ • ■) to name 

t{s) , ^'^'"^ , , . (18) 

Those partial waves with (U) = (00), (11), (20) are used to construct |T|^. Agreement with experimental 
data on pion scattering is very good and the theoretical constraints of unitarity and chiral symmetry are 
properly implemented. 

The low-energy pion resonances in the scalar (a) and vector {p) channels naturally appear in this descrip- 
tion, although not as elementary degrees of freedom, but as dynamically generated resonances. 



6. Results and Discussion 



Our numerical computation for the bulk viscosity of the pion gas, based on elastic processes, is plotted 
in Figure 1 (left panel). We observe a broad maximum at low temperatures, consistent with ^ = at 
zero temperature (a non-relativistic prediction for the monoatomic gas), and a clear growth of ( at larger 
temperatures approaching the phase transition to the quark-gluon plasma. Starting at a higher, 400 MeV 
temperature, we have plotted also the asymptotic estimate for ^ in that phase [5] . In the right panel we also 
give these results quotiented by the entropy density s, since this is a dimensionless number that characterizes 
the sound wave attenuation. For SU(3) and Nf = 3 we have [9], in the high temperature regime 

^=2-10-4 , g-^(T) = -^log(^—]+A^log(2log(^—]]. (19) 

s iog(MM^ ^ ^ ^ 87r2 ^\3QMeVj 9tt^ ^\ ^\30MeYjJ ^ ' 

This ratio we find ourselves in the pion gas to be monotonously decreasing as a function of temperature. 
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Fig. 1. Bulk viscosity of a pion gas in tiic Inverse Amplitude Method (left) and ratio of bulk viscosity to entropy density (right). 
In addition we plot the asymptotic quark-gluon plasma estimate of [5]. 



Since we have contributed for the first time an evaluation that ahows for the chemical potential, we also 
repeat the plot in Figure 1 for various values of /x in Figure 2 (for the pion gas only). 

We should like to remark again that the natural cnd-valuc is fj, = m,r which corresponds to chemical 
equilibrium. Increasing // generally yields a larger viscosity. 
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Fig. 2. Same as in Figure 1 as a function of the chemical potential for the pion gas. 



Since this calculation employs two quantities that are not widely used in the literature (the sound speed 
at fixed particle density as appropriate for the fixed-chemical potential calculation, and the inverse com- 
pressibility), we plot them also in Figure 3 for the reader's convenience. 



We also wish to compare our numerical results with prior approaches based on the elastic pion-pion 
interactions. We have chosen two works. The first [2] employs a pion scattering amplitude that fits the 
experimental phase-shifts but has no connection to chiral perturbation theory. The author does not plot 
the viscosity for low temperatures. Where he provides his data, our calculation is numerically similar but 
somewhat higher. 

The second computation [1] is a field theory evaluation based on a certain ladder resummation, and is 
numerically off our result based on the physical phase shifts. However, the qualitative features and saliently 
the low-temperature limit coincide with our findings. 

Those authors also mention that the first peak is due to the scale provided by the pion mass, m,r- Therefore 
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Fig. 3. Sound speed v„ at fixed particle density (left figure) instead of fixed entropy as conventionally used when n = 0. Inverse 
compressibility (right figure) Kp^ in units of the pion mass. 
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Fig. 4. Comparison of our computation with prior evaluations at /i = 0. The solid line is our computation with the Inverse Am- 
plitude Method at zero chemical potential. The numerical result is similar to Davesne's evaluation based on a phenomenological 
parametrization of the physical pion phase-shifts [2] (dotted line at the bottom). Extending the computation to lower temper- 
atures, we find a maximum and then a return to small values of the bulk viscosity, in agreement with Fernandez-Fraile and 
Gomez Nicola in chiral perturbation theory (that is however numerically a little too large) , and the non-relativistic vanishing 
behavior of C for a monoatomic gas. 

we vary this mass in our computer program, and plot the result in Figure 5. While true that the position 
of the maximum reacts to the change in the pion mass, it does so in an unpredictcd way. Upon increase of 
the pion mass, the bulk viscosity actually decreases, feature not trivially in accordance to the expectation 
of proportionality to the scale-invariance anomaly C oc The position of the peak does shift to higher 
temperatures (qualitatively in agreement with the scale-invariancc violations appearing at T ~ m.^). 
Upon increasing the chemical potential this peak is seen to reduce and even disappear at /i = rUn. We 
temptatively conclude that this peak might be caused by a mismatch and interplay of the two parameters 
/U and m,r- 

To summarize, we have provided a wrap-up of what is known from the bulk viscosity in the low-energy 
pion gas, and find the transport coefficient to be dominated by elastic processes. 
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Fig. 5. Pion mass dependence of the bulk viscosity. The peaJc at moderate temperature is seen to be very sensitive to the pion 
mass. 

With an independent Boltzmann calculation we have extended the computation to include finite chemical 
potential. We find qualitative agreement with prior estimates based on those elastic processes. We look 
forward to heavy-ion experiments being able to access this interesting transport coefficient proportional to 
the trace anomaly in the underlying quantum field theory. In a separate publication [4] we are proposing an 
observable that can be used experimentally. 

Note added: In the final phases of composition of this work we have been made aware of a computation [13] 
of the contribution of inelastic pion scattering processes to the scattering rate for a pion in the medium. 
This contribution to the rate, that does not suffer from the particle-conservation zero mode, is found to be 
exponentially suppressed in agreement with our argument. 

Work supported by grants FPA 2008-00592, FIS2008-01323, FPA2007-29115-E, PR34-1856-BSCH, UCM- 
BSCH GR58/08 910309, PR34/07-15875 (Spain). JMT-R is a recipient of an FPU scholarship. 
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